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vo', 

We develop a chiral symmetric relativistic mean field model with logarithmic sigma 
potential derived in the strong coupling limit of the lattice QCD. We find that both of the 
nuclear matter and finite nuclei are well described in the present model. The normal vacuum 
is found to have global stability at zero and finite baryon densities, and an equation of state 
with moderate stiffness (K ~ 280 MeV) is obtained. The binding energies and charge radii 
of Z closed even-even nuclei are well reproduced in a wide mass range from C to Pb isotopes, 
except for the underestimates of binding energies in several jj closed nuclei. 
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§1. Introduction 



> 

Chiral symmetry is a fundamental symmetry of QCD at zero quark masses, and 
its spontaneous symmetry breaking generates constituent quark and thus hadron 
masses. The Nambu-Goldstone boson of the chiral symmetry, i.e. pion, mediates 
the long range part of nuclear force, ^ and the midrange attractive nuclear force can 
be described with a light scalar isoscalar meson, so called a, which would be the 
chiral partner of the pion representing the fluctuation of the chiral condensate. Thus 
it would be desirable to respect the chiral symmetry in theories of quark, hadron, and 
nuclear physics. Actually, many models and theories of quarks and hadrons such as 
the Nambu-Jona-Lasino model, 2 -'' 3 ' the sigma model, 4 ) and the chiral perturbation 
theories 5 ^ have been constructed on the basis of the chiral symmetry. 

In nuclear many-body problems, relativistic mean field (RMF) models have been 
developed to describe properties of nuclear matter and finite nuclei. The first RMF 
model proposed by Serot and Walecka 6 ) involves the coupling of nucleons to the 
isoscalar vector field, the omega meson, in addition to the coupling to the sigma me- 
son. Later on, non-linear self-coupling terms of sigma 7 ) -10 ) and omega mesons 10 )' 11 ) 
are introduced to obtain better descriptions of nuclear matter and finite nuclei. Ob- 
tained Lagrangians contain a 2 , a 3 , a 4 terms, which remind us of the form of the 
chiral linear sigma model. RMF models have been successfully applied to various 
nuclear many-body problems; the nuclear matter saturation, 6 ) single particle levels 
in finite nuclei including the spin-orbit splittings, 6 ) nuclear binding energies 8 ) -10 ) and 
nucleon- nucleus scattering cross sections and spin observables . 12 ) Furthermore, the 
equations of state (EOS) constructed with these models have been successfully ap- 
plied to compact astrophysical objects such as neutron stars and supernovae. 13 ) -15 ) 
Having these successes and the Lagrangian forms in mind, it would be natural to 
expect that RMF is not only a phenomenological model parameterizing nuclear en- 
ergy functionals, but also a starting point of finite baryon density hadronic models 
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provided that the chiral symmetry is respected. 

Contrary to the expectations above, simple chiral RMF models fail to describe 
the nuclei as a many-body system of nucleons. Since the valence nucleon Fermi inte- 
gral prefers smaller chiral condensate, the normal vacuum jumps to a chiral restored 
abnormal one at a low density below the normal nuclear density in a mean field 
treatment of the linear sigma model. 16 )- 18 ) Thus we need to have a more repulsive 
potential than in the linear sigma model at small values of a at vacuum or at finite 
densities. There are several prescriptions proposed so far to generate this repulsive 
potential in order to cure the abnormal vacuum problem at low densities. One of the 
prescriptions is to include fermion vacuum fluctuations. 17 )' 19 )' 20 ) From the baryon 
one-loop renormalization, non-linear and non-analytic sigma potential terms appear 
and stabilize the normal vacuum. 17 )' 19 ) With quark loops, one obtains similar terms 
as in the baryon loops, and they break the chiral symmetry dynamically. 3 ) When 
one tunes the coefficients in the potential terms generated by the quark loops, it 
is also possible to prevent the vacuum from falling into the abnormal vacuum at a 
small density. 20 ) In both of the simple interaction term from fermion loops 

of the form — <r 4 log<7 2 , where a oc — (qq) represents the light scalar field, gives rise 
to instability at large chiral condensates, then the normal vacuum is not a real en- 
ergy minimum state. This may not be a practical problem in discussing uniform 
nuclear matter if we regard the extrema as the normal vacuum, but in finite nuclei 
the solution of the scalar meson tends to be anomalous when the barrier height is 
low. In addition, we should also include meson loop effects in the baryon loop case, 
and meson loops tend to cancel the non-linear sigma terms coming from the baryon 
loops 21 ) and the normal vacuum would become unstable again at low densities. The 
second approach is to introduce the coupling of sigma and omega mesons. 18 )' 22 ) Due 
to the reduction of the omega meson mass according to the partial chiral symmetry 
restoration at finite densities, the repulsive vector interaction becomes strong then 
the normal vacuum is kept to be stable even at high densities. Quantitatively, how- 
ever, the a 2 uj 2 coupling gives rise to large repulsive vector interaction and the nuclear 
matter incompressibility is found to be much larger than the empirical values. 18 )' 22 ) 
In order to soften EOS in this scenario, we need to include higher order terms such 
as (J 6 and a 8 . 23 ) In Ref. 23), it is demonstrated that we can actually construct an 
RMF Lagrangian which gives a soft EOS, but the coefficients of the higher order 
potential terms are negative leading to the instability at large chiral condensates. 
The fourth way to overcome the abnormal vacuum problem is to incorporate the 
glueball field which simulates the scale anomaly in QCD. 24 )~ 29 ) While it may not be 
justified to include the unobserved glueballs in hadronic models, these models give 
nuclear matter incompressibility in an acceptable range, and they roughly explain 
the bulk properties of nuclei. 25 ) From the symmetry requirement, the glueball is 
conjectured to couple with a in the form of — x 4 log a 2 , where x denotes the glue- 
ball field. This potential term ensures the stability of the normal vacuum in the 
mean field approximation: The mass of the glueball is assumed to be heavy, 1 — 2 
GeV, then it would be reasonable to consider that the glueball expectation value is 
constant (frozen glueball model) in low energy phenomena. Under this assumption, 
the above potential term is divergent at a — > to prevent the normal vacuum from 
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collapsing to the chirally restored abnormal one. 

These approaches are based on QCD inspired effective hadronic or quark models, 
and it is preferable to obtain the chiral potential (energy density as a function of a) 
directly from QCD. At present, it is not yet possible to obtain the chiral potential 
in Monte-Carlo simulations of the lattice QCD, since the quark loop contribution is 
very strong in the chiral (massless) limit. One of the promising directions would be to 
invoke the strong coupling limit of the lattice QCD, where it is possible to perform 
some analytic evaluations of the chiral potential in the chiral limit. Actually, a 
logarithmic a potential term similar to that in the glueball model had been already 
derived in the strong coupling limit of the lattice QCD. 30 ) In the strong coupling 
limit, ^/g'qcD ~^ 0> pure gluonic action disappears and we can perform the one- 
link integral analytically 30 ^ 34 ^ After the 1/d expansion, bosonization, and fermion 
integral, we can obtain the effective free energy with a logarithmic term, — logo -2 . 30 ) 

In this paper, we study nuclear matter and finite nuclei in a flavor SU(2) chirally 
symmetric relativistic mean field model containing a logarithmic potential term of a, 
— logo -2 , derived from the strong coupling limit of the lattice QCD 30 )~ 32 ) at vacuum. 
A phenomenological uj self-interaction term, {uj^uj^) 2 , is also included in the effective 
Lagrangian. Requiring that the pion and nucleon masses and the pion decay constant 
are given, we have four free parameters, m^^g^^gp and c w [a mass, coupling constants 
of nucleons with oj and p mesons, and the coefficient of the (oj^uj 11 ) 2 term). Two of 
these parameters (g^ and cj) are determined to fit the nuclear matter saturation 
point, (po,E/A) = (0.145fm -3 , — 16.3MeV), and other parameters (m a and g p ) are 
determined to reproduce the binding energies of Sn and Pb isotopes. By choosing 
these parameters appropriately, we find that the obtained EOS is as soft as those 
in phenomenologically successful RMF models such as NL and TM models. 8 )~ 10 ) 
Bulk properties (binding energies and charge radii) of proton (sub-)closed even-even 
nuclei are also well explained in a comparable precision to NL1, 8 ) NL3, 9 ) and TM 10 ) 
models. 

There are several works including the coupling of nucleons with the negative 
parity baryons 20 )' 35 ) and works based on the other chiral partner assignment of 
pions (vector manifestations). 36 ^ While these are in promising directions to solve 
the difficulty in chiral RMF models, we stick to a naive assignment of nucleon and 
pion chiral partners in this work. Finite temperature and finite chemical potential 
treatments in the strong coupling limit of the lattice QCD would be also necessary 
to describe the chiral phase transition at high temperatures and densities, but these 
are beyond the scope of this paper. In order to discuss the properties of nuclear 
matter and finite nuclei around their ground states, we expect that the dynamics in 
hadronic degrees of freedom dominates. 

This paper is organized as follows. In Sec.|2J we briefly explain the derivation of 
the a potential term in the strong coupling limit of lattice QCD, and describe our 
effective hadronic Lagrangian. In Sec. [3 we discuss the properties of nuclear matter 
and finite nuclei calculated in the present model in comparison with the results in 
some other models. We summarize our work in Sec. UJ 
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§2. Chiral symmetric sigma potential 



2.1. Logarithmic chiral potential in the strong coupling limit of the lattice QCD 

Chiral potential V a (energy density as a function of a at vacuum) is one of 
the most important ingredients in chiral models; it dynamically breaks the chiral 
symmetry, and the finite expectation value of sigma generate the constituent quark 
and hadron masses. In this paper, we utilize the chiral potential derived in the strong 
coupling limit (SCL) of lattice QCD. 30 )~ 32 ) Here we briefly summarize how to derive 
the chiral potential. 

In the strong coupling limit {g — ► oo), we can ignore the pure gluonic part of 
the lattice QCD action which is proportional to 1/g 2 , and we keep only those terms 
including fermions Sf written as, 

Sf[x,X,U] = x(x)Up(x)x(x + p,)-x(x + p)Ul(x)xix) , (2-1) 

X,fl 

where n^(x) = (— i^o+^H hz^-i^ an( j we eX p ress the action in the lattice unit. We 
consider two species of staggered fermions simulating u and d quarks. After integrat- 
ing out the link variable Uu in the leading order of 1/d expansion and introducing 
the auxiliary fields o~ a s, we obtain the following partition function 



Z = J T>[ X ,x,U] exv(-S F [ X ,X,U]-Su[U]) 
^ J V[x,x] e- 



M af 3{x)V M (x, y)M(y)/3 c 



x,y,a,0 



J ^[x,X,o-]exp(-S a [x,X,o-]) , 



(2-2) 



s<r[x,x,<A= Y a ( x )aP v M 1 ( x ^y) c7 (y)^ + 2 a ( x )apM(x)^ a ■ (2-3) 

Mesonic composites are defined as M. a p{x) = Xa( x )Xp(x), where the superscript 
a denotes color and the subscripts a and j3 show the flavor. The lattice mesonic 
inverse propagator V M (x, y) is given as V M (x, y) = £^ (<W+A + <W~a) / 8N ^- From 
the first to the second line in Eq. (|2-2|) . we have used the one-link integral formula, 
J dUU a bUl d = 5 a dd~bc/N c . The auxiliary fields are related to the expectation values 
of the mesonic composites as {a a p{x)) = — (Vm (x, y)-M. a p(y)}- 

Now we consider static and uniform scalar a and pseudoscalar tv condensates, 
and we substitute the auxiliary fields by the mean field ansatz, 



( \ _ 1 72^ + i£ ( x ) n °) H x )< 



(2-4) 



where e(x) = (— l^o+asi+sa+a^ Since fermions are decoupled in each space-time 
point, we can easily perform the fermion integral. The effective free energy is ob- 
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tained as, 

V x {a,7v) = (tr [o-V^a]} - N c log det a al3 = (V^tx [jltf+Af] - N c log det a a(i 

= N c 4? - N c log 4? {4? = a 2 + Tr 2 ) , (2-5) 

where (• • •) denotes the space-time average, and M represents the meson matrix in 
which the e{x) factor is omitted from o~ a p, 

M=[^\ ! c . o . (2-6) 

While the coefficients of (f> 2 and — log (j) 2 are given in Eq. (j2-5[) , these coefficients 
are fragile. In addition to the freedom in the choice of the lattice spacing and the 
scaling factor for a and 7r to be in the canonical form, it has been shown that the 
baryonic composite contribution modifies the coefficient of 2 . 32 ) Furthermore, two 
species of staggered fermions corresponds to Nf = 8, and the coefficient modification 
may not be trivial when we take Nf = 2. Thus we regard them as parameters to 
obtain physical masses of a and tt mesons, and adopt the following coefficients in 
the later discussions. 

Va° L = V x {a, tt) - c a a = b a 4> 2 - a CT log cf 2 - c a a , (2-7) 
a CT = ^f( m l ~ m l) > b a = ^(ml + ml) , c a = f^m 2 . (2-8) 

After requiring that V a has a minimum at a = f n and fitting the pion mass m n , one 
parameter m a is left as a free parameter. 

Because of the singularity of V a at a — >• which comes from the logarithmic 
term, chiral symmetry restoration is suppressed with this chiral potential. One can 
doubt that this singularity may come from an artifact of the strong coupling limit. 
Indeed, in a finite temperature treatment of the strong coupling limit, where the 
anti-periodic boundary condition of fermions are taken care of, we do not have a 
divergent behavior at a — > 0, 33 ) but we still have a finite negative derivative at a = 
even in the chiral limit when T = 0. This finite negative derivative at a = is enough 
to suppress the full chiral restoration at finite density, because the nucleon Fermi 
integral contribution behaves as p B o~ 2 and we always have a minimum at a finite a 
value. Therefore, we consider that the present chiral potential V a would be a good 
starting point to describe cold nuclear matter and nuclei. For finite temperatures, 
the singularity disappears also in the derivative, and the full chiral restoration will 
take place. In that case, we have to take account of the finite temperature effects in 

V 31 )> 32 ) 

v (7 • 

2.2. Comparison with other models 

There is a variety of possibilities for the choice of the chiral potential V a ((j) 2 ). The 
simplest one is that in the chiral linear a model ((j) 4 theory) proposed by Gell-Mann 
and Levy. 4 ) 

V^ = \^-flf + \m 2 ^ 2 -Umla, A = . (2-9) 
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Fig. 1. Energy density at vacuum. Energy densities at vacuum as a function of a calculated in the 
SCL model (solid curve) are compared with those in the linear a (4> 4 , dotted curve), NJL 2 - 1 ' 3 -' 
(open squares), baryon loop (BL, dashed lines) 19 - 1 , Sahu-Ohnishi (SO, dot-dashed lines) 23 \ and 
TM1 (filled circles) 10) models. 



With this chiral potential, we have only one free parameter m a and the theory is 
renormalizable. The double-well structure in the linear a model can be understood as 
a result of quark-antiquark condensation dynamically generated from the chiral sym- 
metric four Fermi interaction of quarks in the Nambu-Jona-Lasino(NJL) model. 2 )' 3 ) 
The chiral potential in NJL is found to be 



TyNJL 
v a 



m 



N c N f 
4vr 2 



Ga\ 



1 + 



X 



log 



1 + y/l + X 2 



(2-10) 
(241) 



where G represents the coupling of the quark to a and tt mesons. 3 ) In Ref. 3), 
parameter values are fixed to fit f n and m n , resulting in the constituent quark mass 
M q = Gf n = 335 MeV and cut off A = 631 MeV. 

Nucleon one-loop renormalization of the chiral linear a model leads to additional 
potential terms at vacuum. In the chiral limit (m^ = 0), the following interaction 



appears, 



19) 



i 

4^2 



log x s 



i 2 ^ 4 
— + x — -X* 
4 4 



(2-12) 
(2-13) 



In Ref. 19), the value of sigma mass is taken to be m a = 572.8 MeV. 
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In Ref. 23), Sahu and Ohnishi proposed to add higher order terms to obtain 
reasonable incompressibility within the Boguta scenario, 18 ' and the chiral potential 
is given in the chiral limit as follows, 



2 



V «° = ^ 2 ~ + fifsoWU) , (2-14) 

fso(x) = ^ (x 2 -l) 3 + ^(. 2 -l) 4 , (2-15) 
6 8 

with parameters m a = 762.3 MeV, C 6 = -74.4 and C 8 = -2.2. 

For comparison, we also refer here the potential in a non-chiral model, TM1. 10 ) 

KT>) = ^ 2 + ^V + f/, (2-16) 

where ip stands for the deviation from the vacuum, ip = f n — a. Parameters in the 
TM1 model are given as m a = 511.198 MeV, g 3 = 15.3383 and g 4 = 0.6183. 

In Fig. ^ we compare the chiral potential in the present SCL model, V^ CL in 
Eq. (|2-7j) . with those in other models. The chiral potential in the NJL model (open 
squares) describes that in the linear a model (0 4 , dotted line) in the region a < f n , 
giving the foundation of the </> 4 model from quark degrees of freedom. As mentioned 
in the introduction, the nucleon Fermi contribution prefers smaller a values at finite 
density, and in order to keep a from collapsing to the abnormal vacuum (a ~ 0) at 
low densities, we have to have more repulsive potential in the region a < f n than in 
the linear a model. While the baryon loop (BL, dashed line) and Sahu-Ohnishi (SO, 
dash-dotted line) models have potentials repulsive enough at small a values, these 
models are found to have instability at large a values. In the present model (SCL, 
solid line), we do not have any instability at any values of a, and we have stronger 
repulsion in the region of a < /„-. It is interesting to find that the SCL model results 
are very close to that in the TM1 (filled circles) except for the diverging behavior at 

(T-»0. 



§3. Nuclear matter and finite nuclei in chiral RMF models 

Relativistic Mean Field (RMF) approach has been developed as an effective 
theory to describe nuclear matter and finite nuclei in a field theoretical treatment 
containing scalar and vector meson couplings to nucleons. In this paper, we consider 
the following chiral symmetric RMF Lagrangian in which nucleons couple with a, 
7r, u and p fields, 

+ X - (d^ad^a + c^tt • d^iv) - V a (a, tt) 

(3-1) 

W^v = d^ u - dyu^ , (3-2) 
Rfw = d^pu ~ %Pn + 9 P Pv x p u . (3-3) 
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Here we have omitted the photon field, which we include for finite nuclear study. 

In this section, we study uniform nuclear matter and finite nuclei with the La- 
grangian in Eq. Q3T[) together with the logarithmic chiral potential V^ CL in Eq. (|2-7jl . 
We search for an appropriate parameter set, containing the a mass (m a ), meson- 
nucleon coupling constants (g w and g p ), and the strength of the u> self-interaction 
(qj), which explains the properties of symmetric nuclear matter and finite nuclei. 

3.1. Nuclear Matter 

First, we study the EOS of uniform symmetric nuclear matter. We assume that 
the meson fields are static and uniform, then the RMF Lagrangian for nuclear matter 
becomes 

+ 1 -ml^ + C ^+ 1 -m 2 p R 2 -V (7 (a), (3-4) 

which includes a, u and p (written as R) mesons. Here we have omitted the Lorentz 
and isospin indices, u> = u>o and R = p^o, for simplicity. 

The energy density in symmetric nuclear matter as a function of a and u> can 
be written as, 

E/V = g N J J^yjp 2 + M* N 2 (a) + g^p B - ^u 2 - + V a (a) , (3-5) 

where g^ = 4 is the spin-isospin degeneracy factor of nucleons, and M^(cr) = g a a 
is the nucleon effective mass. The equilibrium conditions read 

d(E/V) dV a 

~ 9uPs + = , (3-6) 



da da 

9uPb ~ mluj - c^u 6 = , (3-7) 



d(E/V) 2, . „ , .3 



dp M* N 
(2tt) 3 . I p 2 + m *2 {a) 



Ps=9N I ,Zv, , 1^— • (3-8) 

* 2 
N 

In symmetric nuclear matter, we have three relevant parameters, m CT ,g w and 
c w . When we give m a as a free parameter, then other two are determined to fit 
the saturation properties, (po,E/V) = (0.145 fm~ 3 ,— 16.3 MeV). In Fig. |5J we 
show the nuclear matter incompressibility K = Qp^(d 2 (E /V) / dp 2 B ) as a function of 
m a . We find that the incompressibility is smaller than 300 MeV in the mass region 
460 MeV < m a < 540 MeV, which can be regarded as the allowed region. Especially, 
at around the incompressibility minimum (K ~ 279 MeV at m a - ~ 500 MeV), we 
obtain EOS which is as soft as that in the TM1 10 ) model. 

In Fig. [3 we compare the energy density as a function of a at p B = (0 — 5)po 
in several chiral RMF models. In the (j) 4 model (left top panel), the nucleon Fermi 
integral contribution (left bottom panel in Fig. |3J) is stronger than the repulsive 
potential at a < f^, and the vacuum collapses to the abnormal one at a density 
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Imcompressibility vs m G 
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Fig. 2. Symmetric nuclear matter imcomplessibility K. The symmetric nuclear matter incom- 
pressibility K as a function of a mass (m CT ) in the SCL model is shown. For a given m n value, 
the tj-nucleon coupling constant (gj) and the u) self-interaction strength (c^) are determined 
to fit the symmetric nuclear matter saturation point, (po, E/A), and the incompressibility K is 
obtained as a result of fitting. 

below /5o- 16 )~ 18 )' 22 ) In order to avoid this collapsing, Boguta 18 ) replaced the vector 
meson mass term m 2 u; 2 /2 with that of the o~ 2 lu 2 coupling term, 

."Boguta 1 m ul 2 ,2 (n n \ 

£*j = YJ 2 ' ^ ' 

In the case of no ui self- interactions, c w = 0, the above coupling term gives large 
ui values at small a as, uj = f^g^pg/mla 2 , leading to a large repulsive potential 
at finite densities. Because of this repulsion, nuclear matter EOS becomes stiff 
(K > 600 MeV), and a increases again at around p B ~ 0.27 fm -3 , as shown in the 
top right panel of Fig. |21 

In the left panel of Fig. 01 we compare the EOS in the present model with those 
in other chiral RMF models and the TM1 model. We adopt m a = 502.63 MeV as 
a typical value, which gives K = 279.14 MeV. Other parameters are summarized in 
Table Q] We find that the present model (SCL) gives softer EOS than that in the 
linear sigma model in the Boguta scenario (Bog.), 18 )' 22 ) and its incompressibility is 
comparable to those in Sahu-Ohnishi (SO) 23 ) and TM1 10 ) models. The baryon loop 
(BL) model 19 ) gives the softest EOS in the models considered here. However, the BL 
model has the instability at large a values, and the incompressibility is too small, 
K ~ 110 MeV. 

In addition to the similarity in EOS, a and to expectation values in the SCL 
model behave in a similar way to those in the TM1 model at low densities. In the 
right panel of Fig. |1J we show a and to expectation values as a function of p B . We 
find that the a expectation value decreases to around a ~ 60 MeV at p B ~ p in the 
SCL and TM1 models, while the decrease is smaller in the BL and Boguta models. 
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Energy Density at p B =(0-5) p 
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Fig. 3. Energy density at finite densities. Calculated energy densities at vacuum (thick solid 
curves) and p B = (1 — 5)po (thin solid curves, from down to up as the density grows) as a 
function of a are compared in chiral RMF models; the linear a model (<f> 4 , top left panel), the 
Boguta model (Bog., top right panel) 18 '' 22 ', the baryon loop model (BL, middle left panel) 19 ', 
and the SCL model (middle right panel). The nucleon Fermi integral contribution (bottom left 
panel) and the results in the TM1 model (bottom right panel) 10 ' are also shown. The dotted 
curves show the equilibrium point where the energy density becomes the local minimum in the 
range of a < f n . 



The modifications of a and lo from the vacuum values determine the nucleon 
scalar and vector(-isoscalar) potentials, U s (p B ) = —g a ^-K — c) and U v (p B ) = g^co. 
The bulk properties of nuclei such as the binding energies and nuclear radii, which 
TM models have been modeled to explain, would be mainly determined by the 
nuclear matter properties and meson expectation values at around and below pq. 
Thus the above similarity of the meson fields at low densities together with the 
EOS similarity suggests that the bulk properties of finite nuclei would be also well 
described in the SCL model as in TM models. 

3.2. Finite Nuclei 

In describing finite nuclei, it is numerically preferable to represent the Lagrangian 
in the shifted field ip = f n — a and to separate the a mass term from the chiral 
potential V a , 

V a (a) = l -m a p 2 + V v {<p) , V v (<p) = -20^ (^j (3-10) 
/scl(z) =log(l-x)+x + y , (3-11) 
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Fig. 4. Equation of state of symmetric nuclear matter (left panel) and density dependence of meson 
fields (right panel) . We show the binding energy per nucleon in symmetric nuclear matter (left) 
and the a and ui field expectation values (right) as a function of the baryon density. We compare 
the calculated results in SCL (solid curves) with those in the Boguta (Bog., dotted curve) 18 '' 22 \ 
baryon loop (BL, dashed lines) 19) , Sahu-Ohnishi (SO, dot-dashed lines) 23) , and TM1 (filled 
circles) 10 ' models. 



since the boundary condition is given as <p — ► at r — > oo. In addition, it is necessary 
to include the photon field which represents the Coulomb potential. Here we take 
the static and mean-field approximation for boson fields, then RMF Lagrangian can 
be written as follows, 

£* MF = W ~ M *n&) ~ 1°U v (lo, R, A)] 

- l -(yyf- l -mlv 2 -V^) 

+ \ (V.) 2 + \mU + ^ + \ (VRf + \ml& + \ (V ' Af , (3-12) 
M* N ((p) = M N - g a (p , U v (u,R,A)= gu ,uj + g p T 3 R + e^p-A . (3-13) 
The field equations of motion derived from this Lagrangian read, 

[-ia ■ V + 0M* + U v ] V'jv = £ii>N , (3-14) 
(-A + ml)<p = g aPs -^- , (3-15) 

(-A + ml)uj = g^Ps - c^lj 3 , (3-16) 
(-A + m 2 p )R = g p p T , (3-17) 
-AA = epl , (3-18) 

where p s = f% + p™, p B = Pb + P%-, Pt = Pb ~ Pb denote scalar, baryon and isospin 
densities of nucleons, respectively. Total energy is given by the integral of the energy 
density given as, 

E = ^2 2|K|£iKQ - - J {g*<PPs + 9ujup B + g p Rp T + e 2 Ap p B ] dr 

i,K,a 
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where we use the Eq. ([3-14j) - ([3-18j) to calculate second order derivatives of meson 
fields. We solve the self-consistent coupled equations l|3-14[l - (|3-18)) by iteration until 
the convergence of total energy is achieved. In this work, we assume that the nucleus 
which we treat is spherical, then the nucleon wave functions are expanded in spherical 
harmonic basis as follows, 

Ipaium = ( \tJol I Co 5 (3-20) 



E 



[F?Jr\<l>- 

)(|G?(r)| 2 + |Fr(r)| 2 ) , (3-21) 



Ps = Y;(^)(\ G nr)\ 2 -\F?(rr) . (3-22) 

i 

where C, a represents the isospin wave function, proton or neutron, and k = I (— (Z + l)) 
for l=j-l/2 {1= j + 1/2). 

In comparing the calculated results in mean-field models with the experimental 
binding energies and charge radii, we have to take account of several corrections. 
In this work, we consider the center-of-mass (CM) kinetic energy correction on the 
total energy and CM and nucleon size correction on nuclear charge rms radius in the 
same way as that adopted in Ref. 10). The CM kinetic energy is assumed to be 

where Pqm = Si Pi IS the CM momentum. This correction gives an exact result 
when the state is represented by a harmonic-oscillator wave function, and we assume 
that it also applies to the RMF wave functions. The CM correction on the proton 
rms radius is written as 



5(rl) = -2(Rc U ■ R P ) + (J^ M ) 

3h 



2 AM n lo 

2 -^(Rl) + (Rl M ) = + ^ (for light nuclei) 



(for heavy nuclei) , 

(3-24) 



where R p = J2ie p r i/Z is the proton CM position, and (r 2 ,) and (r^) represent the 
proton and matter mean square radii, respectively. We assume again that harmonic- 
oscillator results applies for heavy nuclei. For light nuclei, we evaluate the correction 
in RMF wave functions, and we consider only the direct-term contributions. The 
charge rms radius is obtained by including the finite size effects of protons and 
neutrons, 



(»&) = ( r D + (^ch)p + f (^ch)n . (3-25) 
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Parameters 

m a (MeV) g p 
502.63 13.02 4.40 200 



Constants 



Mjv (MeV) 


U (MeV) 


m* (MeV) 


m u (MeV) 


m p (MeV) 


9<y = Mn/ fir 


938 


93 


138 


783 


770 


10.08 




he (MeV-fm) 


po (fm" 3 ) 


E/A( Po ) (MeV) 


hu (MeV) 


( r ch)p 




197.32705 


0.145 


-16.3 


41 A~ 1/3 


(0.862 fm) 2 


-(0.336 fm) 2 



Table I. Parameters and constants adopted in the present work. Two of the parameters (g^ and 
Cuj) are determined to fit the saturation point of symmetric nuclear matter and others (m CT and 
g p ) are fixed through global fitting of Sn and Pb isotopes' binding energies. 



We evaluate the binding energies and charge rms radii with these corrections, and 
the pairing energy for open-shell nuclei are neglected. 

In describing finite nuclear properties, we have two free parameters, g p and m a , 
which we cannot determine from the symmetric nuclear matter saturation point. We 
have fitted the binding energies of Sn (Z = 50) and Pb (Z = 82) isotopes and have 
fixed the parameter values as g p = 4.40 and m a = 502.63 MeV. Other parameters 
are obtained to fit the symmetric nuclear matter saturation point, and summarized 
in Table QJ By using this parameter set, we have calculated the binding energies and 
charge rms radii of C, O, Si, Ca, Ni, Zr, Sn and Pb isotopes. 

In Table [HJ we show the calculated results of binding energies per nucleon and 
charge rms radii of doubly (sub-s)closed stable nuclei, 12 C, 16 0, 28 Si, 40 Ca, 48 Ca, 
58 Ni, 90 Zr, 116 Sn, 206 Pb and 208 Pb, in comparison with the experimental data and 
the results in TM, 10 ) NL1 8 ) and NL3 9 - 1 models. We find good overall agreement of 
the SCL results with the experimental data for heavy nuclei. For light nuclei, we 
underestimate the binding energies of nuclei, 12 C, 28 Si and 58 Ni, by 0.2—0.6 MeV. 
These nuclei have proton or neutron numbers of Z{oi N) =6, 14 and 28, i.e., they 
are jj closed nuclei. These underestimates imply that the spin-orbit interaction in 
the SCL model is not enough to explain the splitting. There are still discussions on 
the strength of the spin-orbit interactions in RMF, 37 ) and it is recently suggested 
that the explicit role of pions has large effects in jj closed nuclei. 38 )~ 40 ) Thus the 
underestimate may be related to the explicit pion effects, appeared clearly with the 
restriction from the chiral symmetry. 

In Fig. |S1 we show the binding energies per nucleon in C, O, Si, Ca, Ni, Zr, Sn 
and Pb isotopes. Except for the above mentioned underestimate of the spin-orbit 
splittings and the underestimates in heavy Zr isotopes which would be due to the 
deformation, 10 )' 41 ) binding energies are well explained in one parameter set in the 
present SCL model from C to Pb isotopes. In Fig. [51 we also find that the calculated 
results of binding energies in the SCL model for heavy nuclei are very close to those 
in the TM1 model. Here we compare the SCL model results with those in the 
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E/A (MeV) 



Nucleus 


exp. 


SCL 


TM1 


TM2 


NL1 


NL3 


1/110 


IF/110 


VIIIF/100 


12 C 


7.68 


7.09 




7.68 












ICQ 


7.98 


8.06 




7.92 


7.95 


8.05 


7.35 


7.86 


7.18 


28 Si 


8.45 


8.02 




8.47 


8.25 










40 Ca 


8.55 


8.57 


8.62 


8.48 


8.56 


8.55 


7.96 


8.35 


7.91 


48 Ca 


8.67 


8.62 


8.65 


8.70 


8.60 


8.65 








58 M 


8.73 


8.54 


8.64 




8.70 


8.68 








90 Zr 


8.71 


8.69 


8.71 




8.71 


8.70 








116 Sn 


8.52 


8.51 


8.53 




8.52 


8.51 








196p b 


7.87 


7.87 


7.87 




7.89 










208 p b 


7.87 


7.87 


7.87 




7.89 


7.88 


7.33 


7.54 


7.44 




charge rms radius (fm) 


Nucleus 


exp. 


SCL 


TM1 


TM2 


NL1 


NL3 


1/110 


IF/110 


VIIIF/110 


12 C 


2.46 


2.43 




2.39 












16 o 


2.74 


2.62 




2.67 


2.74 


2.73 


2.64 


2.62 


2.69 


28 Si 


3.09 


3.04 




3.07 


3.03 










40 Ca 


3.45 


3.44 


3.44 


3.50 


3.48 


3.47 


3.41 


3.40 


3.45 


48 Ca 


3.45 


3.46 


3.45 


3.50 


3.44 


3.47 








58 Ni 


3.77 


3.77 


3.76 




3.73 


3.74 








90 Zr 


4.26 


4.27 


4.27 




4.27 


4.29 








116 Sn 


4.63 


4.62 


4.61 




4.61 


4.61 








196p b 




5.48 


5.47 




5.47 










208p b 


5.50 


5.54 


5.53 




5.57 


5.58 


5.49 


5.49 


5.53 



Table II. Experimental and calculated binding energy and charge rms radius values of stable nuclei. 
Calculated results in the SCL model are compared with those in TM1 10) , TM2 10) , NL1 8) , NL3 9) , 
glueball model (I/110) 25) , frozen glueball model (IF/ 110 and VIIIF/100) 25) and experimental 
ones, respectively. 



TM1 model without pairing corrections. As shown in the top panel of Fig. these 
two models give similar results of nuclear densities. The agreement in the density 
distributions is not surprising since the behavior of a and oj mesons as well as the 
EOS is very similar in these two models at low densities, as shown in Fig. |3J In 
the bottom panel of Fig. H3 we find some differences and changes of the order in the 
single particle energies, especially in those for neutrons around the Fermi energy, 
which would be affected by the model details. 

There are several other chirally symmetric models which well describes the nu- 
clear matter as well as finite nuclei of 16 0, 40 Ca and 208 p D- 25 )> 26 ) In these models, 
they introduce the broken scale invariance through the glueball field, and this glue- 
ball couples to a in a logarithmic term, which is conjectured under the requirement of 
the scale invariance. In Ref. 25), Heide et al. proposed a chiral effective Lagrangian 
containing the glueball, and examined that model through nuclear matter and finite 
nuclear studies. With some of their parameter sets (1/110, IF/110), one can describe 
the finite nuclear properties reasonably well, but these parameters give stiff nuclear 
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3 i ^ , , , 1 , , ^ 

10 20 30 40 50 100 150 200 
A A 

Fig. 5. Nuclear binding energies per nucleon of Z (sub-)closed nuclei. Calculated results in the 
SCL model (solid lines) are compared with the experimental data (points) and the calculated 
results in the TM1 10) model (dotted line). 

matter EOS, K > 340 MeV. With the parameter set of VIIIF, one can obtain reason- 
ably soft EOS (K = 267 MeV), but the binding energies of nuclei are underestimated 
by 0.4 - 0.8 MeV per nucleon for 16 0, 40 Ca and 208 Pb. In Ref. 26), Furnstahl et al. 
extended the terms of the chiral effective Lagrangian with glueball having eight free 
parameters. In their work, the binding energies of 16 0, 40 Ca and 208 Pb nuclei are 
well described (up to 0.4 MeV per nucleon deviation from the experimental data), 
and reasonably soft EOS's (K = (194 — 244) MeV) are obtained. Compared to 
these models, the present SCL model has only four free parameters which have been 
determined to fit the symmetric nuclear matter saturation point and finite nuclear 
properties. The main difference in the SCL model from these two glueball mod- 
els is in the uo self-interaction term, c^u; 4 , which simulates the suppression of u at 
large densities, which has been demonstrated to emerge in the Dirac-Briickner HF 
theory 11 ^ 

§4. Summary 

In this work, we have developed a relativistic mean field (RMF) model with a 
chiral symmetric effective potential V a (vacuum energy density as a function of a) 
having a logarithmic term derived in the strong coupling limit (SCL) of the lattice 
QCD. The logarithmic potential term of a is found to have favorable features; it 
prevents the normal vacuum from collapsing at low densities, and it does not have 
any instabilities as a function of a. By introducing vector mesons (u and p) and a 
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CJj 

u 
S 



SCL(n) 

SCL(p) 
SCL(tot) 

TMl(n) □ 

TMl(p) o 
TMl(tot.) 




Exp. TM1 SCL 
P 



Exp 



TM1 SCL 
n 



Fig. 6. Comparison of nuclear density (top) and single-particle energies in the SCL and TM1 10) 
models for 208 Pb nucleus (bottom). Calculated results in SCL model are compared with those 
in TM1 model 10 ' and experimental data in Ref. 42). 



non- linear vector meson self- interaction term ((u^u;^) 2 ) in a phenomenological way 
to fit the nuclear matter saturation point and the binding energies of finite nuclei, we 
have demonstrated that both of symmetric nuclear matter equation of state (EOS) 
and finite nuclear properties are well described in one parameter set containing four 
free parameters. The obtained EOS is comparable to those in phenomenologically 
successful but non-chiral models such as TM1, 10 ) and to that in a chiral model 
containing phenomenologically added higher order terms of o\ 23 ) Binding energies 
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of finite nuclei are also well reproduced in a wide mass range from C to Pb isotopes, 
while the binding energies of several jj closed nuclei are underestimated, suggesting 
smaller spin-orbit interactions in the present model. 

Comparisons of the present SCL model are made with other chiral RMF models 
at vacuum and finite densities in symmetric nuclear matter. We have demonstrated 
that other chiral RMF models have some problems in simultaneous description both 
of nuclear matter and finite nuclei: The baryon loop (BL) 19 ) and Sahu-Ohnishi 
(SO) 23 ) models have instability at large a values, and the linear a model with the 
a 2 uj 2 coupling (Boguta model) gives too stiff EOS of symmetric nuclear matter. 
While the chiral RMF models with glueballs, which represent the broken scale in- 
variance, give reasonably good description both of nuclear matter and finite nuclei, 
introducing unobserved glueballs may not be justified in hadronic models. 

The energy functional in the present SCL model seems to be very similar to that 
in the TM1 model at low densities. This point has been examined in the symmetric 
nuclear matter EOS and finite nuclear properties such as binding energies, charge 
rms radii, nuclear densities and single particle levels. At high densities, the SCL 
model gives a little softer EOS than that in TM1. This softness has some effect 
on neutron star properties and supernova explosion energies, which will be reported 
elsewhere. 

There are several directions which should be investigated further. First, the 
binding energy underestimate problem in jj closed nuclei should be studied in an 
extended framework including explicit role of pions. 38 )' 39 ) Next, it would be desir- 
able to extend the present model to that with the flavor SU(3) chiral symmetry. 
Preliminary works in this direction have been already made, which suggest that we 
can obtain a softer symmetric nuclear EOS (K ~ 210 MeV) due to the coupling to 
hidden strangeness. Thirdly, the finite temperature effects have not been included 
in the effective potential V a adopted in this work. Since the present V a diverges at 
a — > 0, it cannot describe the chiral phase transition at high temperatures. There 
are several works which include finite temperature effects in the strong coupling 
limit of the lattice QCD, 31 )' 32 ) and these works have shown that the effective po- 
tential is smoothed at finite temperatures, and it is possible to describe chiral phase 
transitions. 
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